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Abstract 



We investigate a many-body wave function for particles on a cylin- 
der known as Laughlin's function. It is the power of a Vandermonde 
determinant times a Gaussian. Our main result is: in a many-particle 
limit, at fixed radius, all correlation functions have a unique limit, 
and the limit state has a non-trivial period in the axial direction. The 
result holds regardless how large the radius is, for fermions as well as 
bosons. In addition, we explain how the algebraic structure used in 
proofs relates to a ground state perturbation series and to quasi-state 
decompositions, and we show that the monomer-dimer function intro- 
duced in an earlier work is an exact, zero energy, ground state of a 
suitable finite range Hamiltonian; this is interesting because of formal 
analogies with some quantum spin chains. 
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1 Introduction 

In this article we study a many-body wave function for particles on long 
cylinders. The wave function is the product of the power of a Vandermonde 
determinant and a Gaussian. It arises as a model wave function in the 
fractional Hall effect |Lau2] and is known as Laughlin's wave function, but 
the wave function, or variants of it, play a role in other areas too, e.g., 
rapidly rotating Bose gases (see |LSj and the references therein) and classical 
Coulomb systems ("jellium" or one-component plasma), see for example jsl 
[F] . The function also resembles expressions studied in random matrix theory. 

A full discussion of the quantum Hall effect background is beyond the 
scope of this article, and we should stress that our result does not seem 
to have direct implications in terms of quantized conductance. There are, 
nevertheless, reasons that might make the result interesting in the quantum 
Hall effect context. The key words are Laughlin's argument, Chern numbers, 
incompressibility, topological order, and the plasma analogy. 

Laughlin's argument for the integer Hall effect [Lau2] used the cylin- 
der geometry and gauge invariance in order to show that Hall conductances 
should have integer quantization. It was later suggested that a ground state 
degeneracy, possibly due to translational symmetry breaking, is required in 
order to reconcile this argument with fractional conductances |TW] . The 
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relation between ground state degeneracy and fractional conductance can in 
fact be made rigorous on a torus, in the Chern number approach [ASY] . This 
leads to the question whether Laughlin's state, on a cylinder, considered as 
an approximate ground state, is degenerate or not: this was the initial moti- 
vation of the present article, see |JLS] for further references. Our main result 
is that at filling factor Laughlin's state exhibits indeed p-fold symmetry 
breaking, at all values of the radius, complementing previous thin cylinder 
results jRHllJLS] . 

This shows, in a way, that the use of Laughlin's function is consistent with 
an important ingredient of the theory of the fractional Hall effect, incompress- 
ibility or the existence of a gap above the ground state(s), see the review [BFj 
and also the discussion in [BESj . Indeed, on a cylinder, at non-integer filling 
factor, incompressibility implies translational symmetry breaking |Koj . Let 
us mention, however, that our symmetry breaking result holds at all values 
of the radius and of the filling factor, regardless of whether or not there is a 
gap. 

The symmetry breaking proven here can be read as another illustration of 
the geometry-dependent degeneracy which is considered a hallmark of "topo- 
logical order" . We do not wish to discuss this notion here, and instead refer 
the reader to |HMj and the references therein for a discussion in combination 
with the fractional Hall effect. 

The plasma analogy refers to the observation that the modulus squared 
of Laughlin's function is proportional to the Boltzmann weight for a clas- 
sical system of point charges moving in a neutralizing background. It was 
originally invoked [Lau2j in order to justify that, at not too low filling fac- 
tors, in a disk geometry, Laughlin's state describes a homogeneous liquid. 
Interestingly, this very same argument, when adapted to the cylinder geom- 
etry, suggests a spatial periodicity. From this point of view, the symmetry 
breaking proven here is closely related to results for one-dimensional jel- 
lium jKul IBH lAM] . jellium on a semi-periodic strip at the "free-fermion 
point" [CFSj and at even-integer couphng (SWKj . and jellium on quasi one- 
dimensional tubes with periodic or Neumann boundary conditions, at arbi- 
trary coupling [A J Jj . 

A curious aspect of the quasi one-dimensional jellium is that it inter- 
polates between one-dimensional jellium, which is known to have no phase 
transition [Ku] . and higher-dimensional jellium, which is expected to undergo 
a Wigner crystallization phase transition (related to the formation of vortices 
in a rotating Bose gas). An interesting open question is, therefore, whether 
at low filling factor Laughlin's state has a phase transition as the radius 
is varied. We know that for every fixed filling factor, when the cylinder radius 
is sufficiently small, the infinite volume correlations depend on the radius in 
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an analytic way, and there is exponential clustering |JLSj . Our results imply 
that if there is a phase transition, it cannot manifest itself in a change of 
the state's spatial periodicity. A possible scenario, instead, could be that 
exponential clustering is replaced with algebraic decay of correlations, which 
in turn is related to the question whether the gap of some toy Hamiltonian 
vanishes at some finite value of the radius as the radius is increased. 

The simplest motivation, perhaps, is to consider Laughlin's function as 
a partially solvable, quantum many-body toy problem. It is well-known, 
indeed, that the wave function is an exact ground state of a suitable inter- 
action which encodes that the wave function has zeros of a given order as 
particles get close [HI IPT[ ITKj . Our problem can also be seen as a problem 
for fermions, or bosons, on a one-dimensional lattice, and the interaction 
takes the form 

H= F{ni-n2)F{ki-k2)cl^cl^Cn^Cn^ (1) 

ki+k2=n-i+n2 

for some suitable rapidly decaying function F. Studying Laughlin's wave 
function amounts to studying the ground state of this simple looking Hamil- 
tonian. 



Our main results are the following: when we let the particle number go 
to infinity, at fixed radius, all correlation functions have a unique limit, and 
the limit state has a non-trivial period along the cylinder axis. At filling 
factor the period is p times the period of the filled lowest Landau level; 
equivalently, p times the period of the Hamiltonian H from Eq. ([1]). We 
also show that the state is clustering, and that bulk correlation functions are 
insensitive to the precise choice of the domain of integration. This is akin to 
the accumulation of excess charge at the boundary of Coulomb systems. 

These results extend previous results for fermions on thin cylinders |JLS] . 
They leave open, however, some questions that were answered positively on 
thin cylinders. Most notably, on thick cylinders, we do not know whether 
the symmetry breaking is already apparent at the level of the one-particle 
density, and we do not know whether the correlation length is finite. This last 
question is, via the exponential clustering theorem |NS] related to another 
open question, namely, whether the Hamiltonian H has a gap above its 
ground state (see, however, jSFL"*" for numerical results at filling factor 1/3). 



The symmetry breaking proven here is also closely related to results 
for classical quasi one-dimensional Coulomb systems |AJJj . From this per- 
spective our result is, on the one hand, a specialization of jAJJ] to 1 + 1- 
dimensional semi-periodic strips; on the other hand, our result is an improve- 
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ment in the sense that we examine the full quantum mechanical state and 
we prove uniqueness (up to shifts) and ergodicity of the hmiting state. 

Our proofs follow jJLS] by exploiting the algebraic structure of the wave 
function. Some of the results in jJLS] were proven using a perturbative ar- 
guments in the thin strip limit. We adapt methods from |AJJ] in order to 
show that part of these results are not perturbative but true for all radius 
values. In particular, the associated renewal process has always finite mean. 

Another goal of this article is to cast the algebraic manipulations and 
expressions of infinite volume correlation functions from jJLS] in a form that 
allows for a more explicit expression of the relationship between quantum 
mechanics and an associated probability problem. We will see that the for- 
malism implicitly used in [JLSj does not stand alone, but is actually similar 
to structures encountered in magnetic itinerant electrons and in spin chains. 

Quickly summarized, jJLSj exploited an algebraic property of powers of 
Vandermonde determinants in order to give a representation both of the 
normalization Cat = | |\E' at | P and the correlation functions. The normalization 
was expressed as a polymer partition function, or more precisely, a sum 
over partitions of some discrete volume into non-overlapping intervals with 
multiplicative weights. 

This representation of the normalization is similar to a formula for the 
canonical partition function for magnetic itinerant electrons given by Aizen- 
man and Lieb, Eq. (11) in |ALj . The canonical partition function was written 
there as a sum over partitions {ui} of the number of electrons A^, with the 
interpretation that 

"as far as the 2;-component S is concerned (...) the system is in 
a superposition of states in which the particles form independent 
'cliques' of size rij." |ALj . Remark (1) after Eq. (11) 

The clique structure carries over to a representation of the correlation func- 
tions: for quantum spin chains, the state can is a weighted superposition of 
(quasi-) states, corresponding each to a partition of spins into "random clus- 
ters of total spin zero" ( [AN] ) . The weights define a probability measure on 
partitions, and a good understanding of this probability measure is a starting 
point for a good understanding of the quantum-mechanical state. 

Something similar can be done for Laughlin's wave function. The physical 
meaning of the analogous clique structure is unclear, but we shall try to 
give some intuition in terms of a ground state perturbation series. Let us 
explain right away why the representation is so useful for the infinite volume 
limit. We will write expectations as convex combinations of linear functionals 
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{quasi- states), indexed by partitions X: 



11*^11 V 

The weights define a probabihty measure Pat on partitions. The decomposi- 
tion allows us to reduce the question of convergence of quantum mechanical 
states to the much simpler problem of convergence of probability measures. 
In the limit — )■ oo, the measure Pat has a limit P, and 

||^^||2 > J dF{X)uJx{a). 

The measure P is a p-periodic renewal process and it is mixing. The infinite 
volume state inherits the periodicity and the clustering from the measure P. 

Finally, we show that the solvable monomer-dimer function introduced 
in |JLS] arises as an exact ground state. At filling factor 1/3, the Hamilto- 
nian is obtained from H by restricting the interaction range to l^i — /C2I < 3, 
1^1 —n2\<3. This does not enter the proofs of our main results in any way, 
but is interesting for two reasons. First, the monomer-dimer Hamiltonian 
might be a good toy model for open questions about gaps and incompressibil- 
ity. Second, the monomer-dimer Hamiltonian, as a sum of non-commuting, 
positive, local operators Yls^s^s, resembles spin chains with nearest neigh- 
bor spin singlet projections. This calls for a better understanding of the 
relation of our monomer-dimer structure and the valence bond structures 
and dimerization encountered in spin chains |AKLTt [AN] . 

The article is organized as follows: Sect. [2] states the main results in 
terms of the continuum wave function; the lattice formulation, many-body 
Hamiltonian and monomer-dimer model are presented in Sect. El The re- 
newal structure of the wave function and the quasi-state decomposition are 
explained in Sect. |H In Sect. |5l we explain how methods from quasi- ID 
Coulomb systems help strengthen perturbative results from |JLS] . Finally, 
Sect. |6] weaves together the formalism from Sect. H] and the bounds from 
Sect. Elto prove our results on correlation functions. 

2 Main results 

Consider A^ particles moving on a cylinder of radius R > 0. Particles have 
complex coordinates z = x + iy. The real part a; is a coordinate along the 
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cylinder axis, and the imaginary part < y < 2nR is an angular coordinate 
around the cylinder. We are interested in the wave function 



^<j<k<N \ k=l 

(2) 

The parameters £ > and p G N are considered fixed. In the quantum 
Hall effect, £ oc 1/ is the magnetic length, and 1/p is the filling factor. 
When p is odd, is a wave function for fermions; when p is even, \E'jv is 
a wave function for bosons. The multiplicative constant is irrelevant for 
correlation functions; a convenient choice is nevertheless given in Eq. (|T2|) . 

Before we state our results, we mention an illuminating alternative form 
of the wave function's modulus. Using the factorization 

exp{zk/R) - exp{zj/R) = (l - exp([zj - Zk]/R)j exp{zk/R). 

and a completion of squares, we find 



|^'^(2;i,...,2;7v)|2 cx Yl |l 

l<j<k<N 



II 11 gjjp| 



2p 




The Gaussian clearly favors a period pi'^/R in the x-direction. The prefactor 
should not destroy the periodicity, since in the sector xi < X2 < ■ ■ ■ < xn, 



Zj - Zk , , ^ f ( \Xj -Xk\\\ 

exp^ = l + 0^exp( 



is almost unity when Zj and Zk are far apart. Moreover, since the Gaussian 
decays exponentially outside a finite cylinder of length pNi"^ /R, the exact 
choice of the domain of integration A should not affect bulk correlation func- 
tions. Theorems 12.11 and 12.21 below state that this picture is correct; that 
this is true for all values of the cylinder radius, and not just for a radius 
small compared to the magnetic length is in essence the main result of this 
paper. 

Remark. The dimensionless parameters are p G N and the ratio 7 := l/R. 
In the proofs, we will choose length units such that £ = 1, so that the period 
of the state becomes pl'^ / R = J97. 
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We now state our results in more detail. Let Z = Wx [0,27ri?]. For a 
given domain of integration A G Z and G N, we consider the ra-point 
correlation function 



Pn ' (^li • • • 5 ^n, ^1; • • • ; ^n) 

N\ 1 



{N-n)\ 11^ 



n\\a Ja 



We wilU consider domains of the form 

A = [a, h] X [0, 2txR], -oo < a < C, Npf/R - C < b < oo (4) 

for some A^-independent constant C, but a and b possibly A^-dependent. For 
example, we can integrate over the infinite cylinder Z : — oo < x < oo, the 
semi- infinite cylinder x > 0, or a finite cylinder of length pNi^/R. 

We are interested in the limit -> oo at fixed cylinder radius R > 0, 
and fixed parameters i, p. The cylinder length 6 — a as in Eq. goes to 
infinity. In this limit, the bulk correlation functions have a unique limit, and 
the limit is p£^/-R-periodic. More exactly, we shall see that 

N N 

(U'^^py,)) Pn{{z,-pf/R}U; {z'^-piVR}]^,) 

= Pn{zi, . . . , Zn', ^1, . . . , Z^). (5) 

for all Zi, . . . , z'^ E Z, and u{py) = exp{ipy /R); the unitary u{y) comes from 
the magnetic translation 

t{f/Re^)tlj^ (z) = exp{iy/R)tP{z - f/R). 

Theorem 2.1 (Infinite volume correlation functions, p-periodicity) . Let N — )■ 

oo at fixed R and A as in Eq. . Then there is a unique family of infinite 
volume correlation functions Pn{{zj}] {zj}) such that, for suitable e{d) with 
e{d) — )■ as d ^ oo, and suitable constant Kn, 

Pn ' {Zlj ■ ■ ■ J Zn', ^1, • • • 5 Z^) — Pn{zi, . . . , Zn', Z^, . . . , Z„) < KnSid), 

whenever ^i, . . . , G A satisfy Xj, x'j > d and Xj, x'j < pNf'/R — d for all j . 
The infinite volume correlations are pi"^ / R-periodic in the sense of Eq. ([5]). 
They do not depend on the precise choice of A. 
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The error terms Kn and e{d) may depend on C in the choice of A, Eq. (jl]). 
The infinite volume correlation functions are invariant (in the usual sense) 
with respect to arbitrary shifts in the |/-direction, yj yj + s, and with 
respect to reversals zj — )■ —Zj. 

Remark. The correlation functions pn determine uniquely a state (■) on a 
suitable bosonic or fermionic observable algebra, and we also have conver- 
gence of states on that algebra, i.e., (\I^7v, a\I^7v)/||^Ar|P — > (a). 

The periodicity statement in the previous theorem does not yet guarantee 
that pi'^/R is the smallest period of the state, but the next theorem does. 

Theorem 2.2 (Symmetry breaking). For all i,R > and p G N, there 
is some correlation function Pni^i, . . . , z^, zi, . . . , Zn) that has p£'^/R as its 
smallest period. 

This should be contrasted with the result from |JLS] : on thin cylinders 
{R i), we know that p£'^/R is the smallest period of the one-particle density 
Pi{z; z), which is a function of x alone. Theorem 12.21 leaves open whether this 
stays true for thick cylinders. In any case, however, the one-particle density 
has a non-trivial period which must be a multiple of the filled lowest Landau 
period i^/R [JLSl Lemma 2]. 

The use of higher order correlation functions in the formulation of sym- 
metry breaking may look confusing, so let us try to make it more concrete. 
Suppose we are given samples of electron gas in the state \I/Ar, and we can 
make repeated measurements of the number of particles N{x), N{x) in small 
width e annuli at two different abscissas x and x. This gives two histograms 
which will look different unless x — x is a multiple of pi'^/R. On thick cylin- 
ders, the histograms might have the same average {N{x)) = {N{x)) but 
different shapes, so we might have to look at higher moments (A^(x)™) in 
order to see the difference. 

The C* algebraic content of Theorem 12.21 is that the states ui, . . . ,Up as- 
sociated with the bulk correlation functions (pn) and its shifts by i'^/R along 
the cylinder axis are distinct. We shall see that they are actually disjoint in 
the sense of [BRl Sect. 4.2.2, p. 370], see Sect. 16.51 This corresponds to the 
mutual singularity of probability measures proven in [A J J] . 

The correlation functions are clustering. We formulate the next theorem 
for integration on the infinite cylinder Z, but note that the insensitivity to 
the precise domain of integration stated in Theorem 12.11 allows us to transfer 
the statement to other domains of integration. 

Theorem 2.3 (Clustering). Let m, n G N and {zj,Zj; 1 < j < m}, {zj,Zj : 
m + 1 < j < n} be coordinate clouds having distance > d > to each other. 
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Then 




' m.+n 



for some N -independent function e{d) with e{d) — as d ^ oo, and some 
suitable constant K^+n- 

On thin cylinders, it was shown in |JLS] that there is exponential clus- 
tering, i.e., e{d) — )■ exponentially fast. Whether this stays true on thick 
cylinders is an open question. 

Clustering is not only of interest in itself, but also enters the proof of 
the insensitivity of bulk correlation functions to the precise choice of domain 
of integration. We will see, indeed, that a switch from integration on the 
infinite cylinder to integration on a finite cylinder can be seen as a perturba- 
tion localized at the boundary; clustering allows us, then, to show that this 
boundary perturbation does not affect bulk correlations. 

Theorem 12.31 is a statement about diagonals: in each of the coordinate 
clouds, the number of primed variables equals the number of unprimed vari- 
ables - we do not separate Zj from Zj. There is off-diagonal decay too, but 
do not offer any explicit estimate, except for the one-particle matrix. 

Proposition 2.4 (Off-diagonal decay of the one-particle matrix). For some 
suitable constant K and all N , z, z' , 



Prop. 12.41 is proven in Sect. 15.31 

3 Interacting fermions and bosons 

For the proofs, it is convenient to view the wave function as a wave function 
for fermions, or bosons, on a one- dimensional lattice. Here we explain how 
this is done. In addition, we recall the expression of the many-body Hamilto- 
nian whose exact ground state Laughlin's function is, and we show that the 
monomer-dimer function from |JLS] is the exact ground state of a suitable 
Hamiltonian too. 

One-dimensional lattice Straightforward algebra shows that \1/Ar(2;i, . . . , zjq) 
is a linear combination of products of one-particle functions 
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with < k < pN — p ; see Eq. f[T^ below. The proportionahty constants are 
chosen positive and such that the ipk{z), /c G Z, are orthonormal in L'^{Z). 

The one-particle function ipk{z) is localized in the cylinder axis direction 
around x = ki'^/R. Thus we may identify lattice sites /c G Z with the orbitals 
ipk{z) and view our wave function as a function for bosons, or fermions, on a 
one- dimensional lattice. 

It is convenient, therefore, to work with creation and annihilation oper- 
ators and Cfc, whose definition we briefly recall. When p is odd, let F be 
the fermionic Fock space associated with the one-particle space L'^{Z). For 
k & Z, the creation operator cl is the linear map in given by c^F = ipk/\F, 
with the antisymmetrized product / A g{z,w) = [f{z)g{w) — g{z)f{w)]/\/2. 
The annihilation operator Ck is the adjoint of in J-". When p is even, 
the definitions are similar, except the wedge product is replaced by a sym- 
metrized tensor product and we have to work in bosonic Fock space. Thus 
creates a particle at lattice site /c, or in the orbital ipk{z), and Ck annihilates a 
particle. The number of particles at site k is given by the occupation number 
operator hk := c^c^. 

We are going to investigate expectations 

(^jv,a^Jv)A 

Wn,A ■■- ||2 • 

II^jvIIa 

for observables a that are products of creation and annihilation operators 
and Ck- The continuum n-point correlation functions is easily recovered from 
expectations of such operators. For example, the one-particle function is 

Pl'^{z;z')= {c*^Ck)N,Ai^k{z)lpm{z') 

In our case, the expectations (c^Ck) n,a actually vanish when k ^ m. 

When p = 1, the wave function is a simple Slater determinant and de- 
scribes independent fermions. As soon as p > 2, it has built-in correlations 
and is the ground state of a suitable many-body Hamiltonian. 



Many-body Hamiltonian The wave function \i/Ar at filling factor 1/p has 
a zero of order p when two particles get close, zj — Zk — )■ 0. In the subspace 
spanned by products of one-particle functions 'ipoi^z), . . . , 'ippN-piz), \i/Ar is the 
unique wave function with this property. This leads to the characterization of 
\E'Ar as the unique ground state of a suitably defined Hamiltonian [H |IPTt[TK] . 
The aim of this section is to recall the form of this Hamiltonian for the 
cylinder problem, see |LH iSFlT^ for similar expressions at filling factor 1/3. 
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This will not be used in any of our later proofs, but helps the interpretation 
of our results in terms of a quantum many-body problem. 

Let Hn{z), n G No, be the Hermite polynomials, given by the generating 
function relation 

exp{2tz -z^) = y2 -,Hn{z) 

n=0 

and Fnit) := exp(-tV4). At filling factor let F (t) := Fo{t) + ■ ■ ■ + 

Fp-i{t)] alternatively, sum only those Fk{t) where < A; < p — 1 and k has 
the same parity as p. For example, F{t) = Fi{t) when p = 3. Consider the 
formal sum over ki, k2,ni,n2 G Z with ki + k2 = rii + n2'- 

H= Yl F{{n,-n2)i)F{{k,-k2)i)cl^cl^Cn,Cn,. (7) 

(Recall 7 = l/R.) We define the finite volume Hamiltonians via free bound- 
ary conditions: For L = {0, . . . ,pN — p} we let He be the Hamiltonian ob- 
tained from this formal sum by restricting summation to integers /ci, /c2, ^i, ^2 G 

Proposition 3.1. The Hamiltonian He is positive He > 0, \&7v ^-5 exact 
ground state Hc^ n = 0, and is the unique ground state of He, considered 
as a Hamiltonian for N particles in the finite lattice C = {0, 1, . . . ,pN — p}. 

The content of the Proposition seems is well-known in the physics lit- 
erature. For mathematical completeness, we include nevertheless a proof 
adapted to the cylinder setting. 

Proof. The positivity of the Hamiltonian follows from an alternative expres- 
sion of Eq. ([7]) 

H= Yl ^*sBs, B,= Y^ Fi2k-i)cs.kCs+k. (8) 

s6(l/2)Z k 

see [LLt Eq. (2)]. The first sum is over half-integers s, the second sum is 
either over integers G Z (when s is integer) or half-integers k G ((1/2)Z)\Z 
(when s is half- integer but not integer); note that the sum of two integers 
and their difference have the same parity, 2s = ki + k2 = ki — k2 mod 2. 
For the ground state property, let us first look at two-particle functions. Let 

^(2^1,^2)= ^ A(mi,m2)V'mi(2;i)^/'m2(^2) 
mi ,m2 
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with the normahzed one-particle functions ipmi^) oc exp(— m^7^/2) exp{im'-fy— 

We use the generating functions relations for the Hermite polynomials: 
and expand 

n=0 

where 

mi+m2 

mi ,r?i2 

Looking at the norm of $n(^) we see that $n(^) = if and only if 

^ A(mi,m2)F4(mi - 7712)7) i^n((A:i - k2h)X{h, = 0. (9) 

mi+m2=ki+k2 

^(^1,2^2) has a zero of order > p as two particles get close if and only if 
^o{Z), . . . , vanish identically, which because of Eq. (Q happens if 

and only if H£^>) = 0. Note that, depending on the parity of p, symmetry 
or antisymmetry gives $2„,(2') = or $2„+i(2') = for free. This proves 
Prop. 13.11 for N = 2 particles. For more than 2 particles, the previous 
argument is applied separately for each pair of particle coordinates {zj,Zk), 
j <k. □ 

Monomer-dimer function In [JLSl Sect. 2] a solvable model was in- 
troduced; this was motivated, roughly, by neglecting the overlap between 
next-nearest Gaussians, or replacing Gaussians with functions of compact 
support. Here we explain an alternative motivation: it turns out that the 
monomer-dimer function is an exact ground state of a Hamiltonian obtained 
from H by keeping only leading order terms. 

We only look at p = 3. We start as in Eq. ((Tj) but keep only terms with 
1^1 — ^2! < 3 and \ni — n2\ < 3. This gives, up to a factor exp(— 7^/2), the 
monomer-dimer Hamiltonian 

jjMD ^ (fikfik+i + 4exp(-37V2)r7fc^fc+2 + 9 exp(-47^)nfcnfc+3 

k 

+ 3exp(-27^)(4+i4+2Cfc+3Cfc + 44+30^+20^+1)). 

Equivalently, 

H^^ = 4exp(-37V2)non2 

+ (4C2 + 3 exp(— 27^)cqC3) (c2Ci + 3 exp(— 27^)0300) + translates. (10) 
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These sums are, of course, formal; given a finite lattice L = {0, 1, . . . , 3A^ — 3} 
we associate as before a Hamiltonian H^^ via free boundary conditions. 

Note that in a thin strip (large 7) limit, truncating the interaction range 
amounts to keeping only leading order contributions, i.e., we neglect high 
powers of exp(— 7^). 

We define a many-body wave function as follows: for /c G Z, let 

%•} := C3fc, ^{fe,fc+i} := -3exp(-272)c^^^^c^^^2- 

Thus a monomer operator creates a single particle and a dimer operator 
creates a pair of particles. The wave function is a sum over monomer-dimer 
partitions (Xi, . . . , Xd) of {0, 1, . . . , — 1} (labeled from left to right): 

"^T ■= Yl I vacuum). 

(Xi,...,Xd) 

"^^^ is the sum of |0 3 6 ■ ■ ■ (3A^ 3)), corresponding to a partition consisting 
of monomers only, and orthogonal terms obtained by hopping one or several 
pairs of particles, for example, 3 6 — )■ 4 5. In particular, all particle pairs have 
distance 1, 3 or higher; there are no pairs of particles at mutual distance 2. 

We have an analogue of Prop. 13.11 We do not, however, prove uniqueness 
of the ground state. 

Proposition 3.2. The Hamiltonian is positive, > 0, and '^^^ is 
an exact, zero energy, ground state. 

Proof. The positivity follows from Eq. flTU]) . In order to check that \l/5v^ ^ 
zero energy state, we first note that has no pairs of particles at mutual 
distance 2 and therefore njnjj^2^^^ = for all j. It remains to see that 

5,*^° = 0, B, := 9+19+2 + 3e-2^'c,c,+3 

for all j. When j is not a multiple of 3, Bj'^^^ vanishes because the wave 
function has no particle pairs at positions (j, j + 3) or (j, j + 1). When j = 3k 
for some k = 0, . . . , N — 1, 

(c3fc+iC3fc+2 + 3e"^'^'c3fcC3fc+3)^Xi " " " |vacuum) 7^ (11) 

implies that either the partition has a dimer {k, k + 1} or it has two monomers 
{k}, {k + 1}. As a consequence we can factorize 

B3k'^T = B3k^{4k4k+3 - 3e"^^'c^fc+ic*fc+2)^|vacuum). 

£ is a sum over partitions of {0, . . . , A; — 1} and 7?. is a sum over partitions of 
{k + 2, . . . , N — 1}. Eq. f lTTj) follows from the observation that i?3fc commutes 
or anticommutes with C and TZ, and 

B3k{4k4k+3 - ^e-^^' 4k+Ak+2) [vacuum) =0. □ 
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4 Renewal structure 



In this section we present the structure of the wave function that is key to 
our proofs. This structure is already visible at the level of normalization 
constants, provided, however, we make a good choice for the multiplicative 
constant k^v in Eq. ([2]). From here on we fix 

1 1 / 1 ^ \ 

(recall 7 = £/R). With this choice, the normalization = W'^nW'z becomes 
supermultiplicative |JLSt Sect. 3.2]. 

In the remaining part of this article we choose length units such that the 
magnetic length is ^ = 1. 



4.1 Block structure of the Vandermonde matrix 

Laughlin's wave function involves the p-th power of a x Vandermonde 
determinant in the variables Zj = exp{'yzj). For 1 < k < N, the upper left 
k X k block of the Vandermonde matrix is itself a Vandermonde matrix. The 
lower right (A^ — k) x [N — k) block is a Vandermonde matrix in which the 
row with variable Zj has been multiplied with Zj-. 



zi 



7k-l 



\ 



ryN-1 

^k+1 



(13) 



Thus the Vandermonde matrix has a simple block structure. Now, the power 
of the Vandermonde determinant is a polynomial and as such can be written 
as a sum of monomials. The block structure suggests that the expansion 
coefficients have some recursive structure. This is, indeed, true, as was shown 
in [FGIL] . Moreover, the recursive structure carries over to coefficients in 
Laughlin's wave function |JLSl Lemma 3]. 

Let us explain this in some more detail. Let ipk{z) be the orthonormal 
basis functions from Eq. We can write 



^7V(^1, 



1 



ttNimi, . . . ,mN)^pmi{zi) ■ ■ -Iprnj^i^ZN)- (14) 



mi , 



,mjv 
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The sum is over integer mi, . . . , m^v running from to pN — p. Not all such 
m's contribute to the sum: if ais[{m) ^ , then for all /c = 1, . . . , A^, 

mi H hmfc>0+pH V p{k - I) (15) 

with equality for k = N. We call such m's admissible. If (mi, . . . ,mAr) is 
admissible and in Eq. (1151) there is equality for some k < N — 1 then 

mi, . . . ,mk < pk - p, pk <mk+i,. . . ,171^ 

and 

aAr(mi, . . . , ttin) = afc(mi, . . . , mk)aN-k{mk+i - pk, . . . ,mN - pk). (16) 

This product rule mirrors the block structure of the Vandermonde matrix f[T51) . 
Not every a]sf{m) factorizes as in Eq. flTBl) : some m's are irreducible. Others, 
in contrast, may factor along more than one k. 

Example. For p = 3, the one- and two-particle wave functions are 

^1 = V'o, = ^0 A - 3 exp(-272)^^i A 7/^2- 

02(1,2) does not factorize, but 02(0,8) does: 

02(0, 3) = 1 = 1 X 1 = oi(0) X oi(3 - 3). 

Remark. Product rules have been derived for other fractional Hall effect trial 
functions [BeR] . 

4.2 Ground state perturbation series 

It would be nice to have a physical interpretation of the product rule de- 
scribed in the previous section. Unfortunately, such an interpretation does 
not seem readily at hand. We hope, however, that the ground state per- 
turbation series argument presented in this section provides some physical 
intuition. The reader should compare the product rule flTB]) to a "ground 
state concatenation" equation, Eq. (fT9|) below. 

In this section we consider only filling factor 1/p = 1/3. We start with 
the Hamiltonian H, or its finite lattice version He, defined in Sect. [31 For 
large 7, the dominant contribution to the Hamiltonian comes from nearest 
and next-nearest neighbor repulsion 

He = exp(-7V2)nfcnfc+i 4 exp{-2-f^)hkhk+2 

kec kec 

+ O(exp(-107V4)). 
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(Thus we discard even more terms than in the definition of the monomer- 
dimer Hamiltonian.) An exact ground state of the nearest and next-nearest 
neighbor repulsion is, obviously, the so-called "Tao-Thouless state" |TT 

= Co C3 ■ ■ ■ C3^„6 4n-3 [vacuum) 

which is why we will abbreviate the corresponding Hamiltonian Hj^. For 
£ = {0, 1, . . . ,pN — p}, \E']7" is actually the unique ground state of Hj^. 
Thus starting from 

Hc^N = {HY + Vc)'^N = 0, HY^Y = 0, 

we obtain 

00 

= ^Y + 5^(-i)"(g(^Z^)'^gv£g)"^]7. (ir) 

n=l 

where Q = Qn,c projects onto the orthogonal complement of "^J/ and 
{HY)~^ is the inverse of HY restricted to TZ{Q). For sufficiently large 
7, the expansion (|T7|) should be convergent. Eq. (fT7|) is a perturbation se- 
ries for the ground state, simplified because both the unperturbed and the 
perturbed ground state energy are equal to 0. 

In the occupation number basis, the Tao-Thouless Hamiltonian HY and 
the projection Q are diagonal, while Vc has non-diagonal, "hopping" , terms, 
e.g., CQC3C2C1. Thus Eq. (fT7|l tells us that the ground state is a sum of the 
Tao-Thouless state and orthogonal terms obtained from "^J/ by hopping a 
pair of particles, (rii, 722) — )■ (^i + r,n2 — r). 

More precisely, let Be be the set of pairs ((/ci, /C2), (ni,n2)) of integers in 
£ with ki + k2 = rii + 77-2, discarding pairs of the form (/ci, ^2) = (^1,71,2) = 
{k, k + 1) or = {k, k + 2) , or permutations thereof. Let 

Acib) := -f{k^ - k2)fin, - n^) QiHYY^Q cl^cl^c^^Cn, Q 

so that 

Q{HYr'QVcQ = -J2 Mb). 

beBc 

Thus \E'Ar is a sum over bond-paths (61, 62, • • • , K) of variable length n, 

00 

^^I^ = ^Y + J2 Yl Acih) ■ ■ ■ Ac{bn)^V ■ (18) 

n=l bi,...,bneBc 

Now let k e {1,...,A^— 1} and imagine drawing a vertical line at a; = pk— 1/2 
splitting the Tao-Thouless state in two. A path (61, ... , 6„) crosses the line 
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if one of the bonds 6j involves integers both to the left and to the right of 
the line. If in Eq. flTSl) we discard crossing paths, the sum factorizes as 



(or a symmetrized product instead of A, depending on the parity oi p). This 
is a /c-particle ground state to the left of the separation line, concatenated 
with an [N — A;)-particle ground state to the right of the line, and is a version 
of the product rule (ITB]) . 

Remark. The idea of treating the full Hamiltonian as a perturbed version of 
a thin cylinder Hamiltonian is related to "one-dimensional" approaches to 
the quantum Hall effect, see [BKj and the references therein. 

4.3 Renewal points and lattice partitions 

Here we explain how the product rule f|T6l) was exploited in jJLSj . First we 
need some notation. We call the points along which ajv(m) factorizes renewal 
points. More precisely, let m = (mi, . . . ,m7v). If m is not increasing, let a 
be a permutation that rearranges the integers, ^^-(i) < ■ ■ ■ < mo-(Ar). We call 
r = pk, k = 0, . . . , N, a. renewal point of m if TT^a{j) = YIiPU ~ '^^^ 
end points and pN are always renewal points. The set of renewal points of 
m is denoted TZ{m). 

Sometimes it is more convenient to work with partitions rather than 
renewal points. Thus let Vn be the set of partitions of the discrete vol- 
ume {0, ...,pN — 1} into discrete intervals or "rods" Xi, of the form 
{pj, . . . , pj + pn — 1}. The number D of rods in the partition varies from 1 
to A^, and the rods are always labeled from left to right, i.e., the elements of 
Xi are smaller than those of X2, etc. We use the short-hand 



and let 71{X) be the set of starting points pj of X^'s, with pN added. We 
write X{m) for the partition associated with the renewal points of m. 

It is useful to group together m's that give rise to the same partition. 
Thus we write 



non-crossing paths 



X — (Xi, . . . , Xd) 




(20) 



with 




X(m)=X 



aAr(mi, . . . ,mN)i^„iAzi) ■ ■ ■ ^JmA^N)- (21) 
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As a consequence of the product rule, the function u;f (zi, . . . , zjy) factorizes 
into an (anti-)symmetrized product of functions Vx associated with individ- 
ual rods instead of the whole partition: if A" = (Xi, . . . , Xd) and p is odd, 

ux{zi, . . . ,zn) = vx^ A ■ ■ ■ A Vxd{zi, . . .,Zn)- (22) 

Good orthogonality properties ensure that 

II*nIIa= Y1 II"'^IIa= Y1 II^^iIIa X ■■■ X ll^xolli 

x&Vn (Xi,...,Xo)e-Pjv 

whenever A has the form [a,b] x [0,27ri?]. The family of functions (vx) has 
the following properties when p is odd: 

• Vx is an antisymmetric function of N{X) = n complex variables, with 
pn the cardinality of X, i.e., X = {pj, . . . , pj + pn — 1} for some j. 

• vx{zi, . . . , zx(x)) depends on X, p and 7, but not on N. 

• Vx is a. linear combination of antisymmetrized products of the ip^s with 
localized indices k: 

Vx e A^(^)span{^fe \ k e X}. (23) 

• Up to a y-dependent phase, vx+p{{zj}) is the same as vx{{zj — P7}). 
Put differently, the shift of a rod is equivalent to the magnetic transla- 
tion of the associated function. (Here p+X could be, e.g., 3+{0, 1,2} = 
{3,4,5}.) 

Similar statements hold for even p, provided we replace antisymmetrized 
products with symmetrized products, e.g., in Eq. (l22l) . 

For later purpose we note that the translational covariance allows us 
to define functions vx for intervals X = {pj,pj + I, . . . ,pj + pn — 1} not 
necessarily contained in {0, 1, . . . ,pN—l}, and functions Ux{zi, . . . , z^) when 
A" is a partition of some set {pj,pj + I, . . . ,pj + pN — 1}, j G Z. 

The decomposition of at as a sum over partitions was exploited in [JLSj 
in order to deduce statements about correlation functions, or expectation 
values of observables. We will essentially follow this approach. First, how- 
ever, let us explain how to go from a decomposition of the vector to a 
decomposition of the state \'^x){'^n\- 
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4.4 Quasi-state decomposition 

To each partition X = (Xi, . . . , Xd) assign the weight 

= Tufr^\M\^ X ■■■ X \\vxn\\^. (24) 

Because of Eq. (14.31) the weights sum up to 1 and define a probabihty measure 
Pat on Vn- It is natural to try to write the state |\l/Ar) {^n\ as a weighted sum 
of states ux, associated each with a partition. Eq. (!20|) immediately yields 

|*jv)(^iv| = Yl MM- 

The state |\l/Ar)(\l/Ar| is, therefore, a sum over pairs of partitions. In order to 
reduce this to a sum over single partitions, we group pairs according to their 
common renewal points. For a given partition X G Vn, let M.i\f{X) be the 
set of pairs {y, Z) such that 

n{y)r\n{z) = n{x). 

Thus in particular, {X, X) G A^Ar(A'). Let ujx be the operator 

" " {y,z)&MN{x) 

By a slight abuse of language we use the same letter for the linear functional 
on operators, 

uJx{,a) = {uz,auy). 

\\Ux\\ 

The sum contains, on top of the diagonal term \ux){ux\, off-diagonal terms 
where y ^ Z. We can write, then, 

(■)7V = = E Vn{?C)ujx. (25) 

which is the decomposition we had looked for. 

Example. Consider the two-particle wave function at filling factor 1/3. Thus 
N = 2 and p = 3. The discrete volume {0,1,..., 5} has two admissible 
partitions: V2 = {V, W} with 

V= ({0,1,2}, {3,4,5}), >V=({0,1,...,5}). 
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We have '^2 = uy + uy^ with 

Mv = ^oA^3, Mw = -3exp(-27^)^i A ?/;2. 

The corresponding probabihty distribution is 

(V) = 1 ^ 9exp(-47^) 

l + 9exp(-472)' ^2(KVj- i^g^^^^_^^2y 

The state decomposes as (■)2 = P2(V)a;v +P2(VV)c<;vv; with 
CUV = |03)(03|, 

cuw = |12)(12| - iexp(272)(|12)(03| + |03)(12|) (26) 

where we have used the short-hand |03) = \ipo A ip^). 

An important remark, immediate from Eq. fl2B]) . is that the matrices ux 
are in general not density matrices, and in general we do not have | Icj^f (a) 1 1 < 
||a||. The picture changes, however, if we evaluate only diagonal operators: 
let a be a product of occupation numbers fik, or any other operator which is 
diagonal in the basis of Slater determinants ipmi A ■ ■ ■ A ifjmN ■ The functions 
ux have a strong orthogonality property which ensures that 

(^a diagonal and y 7^ =^ {uy, auz) = 0. 

As a consequence, if a is diagonal, 

uJx{a) := tT{uJxa) = {ux,aux)/\\ux\\^- 

Thus as far as diagonal operators are concerned, ux is a. state. In the lan- 
guage of [ANj . ojx is a quasi-state: a linear map on a C* algebra whose 
restriction to a commutative sub-algebra is a state, and Eq. fl25l) is a quasi- 
state decomposition. 

The quasi-state decomposition allows us to shift our focus from the com- 
plicated quantum mechanical state to the probability measure Pat, which is 
a much simpler object. This reduction is useful because the ux themselves 
are reasonably simple, as ensured by the following properties: 

• Clustering: let a, b be local observables. "Local" refers to the lattice 
picture: a and b are, for example. Wick-ordered monomials, and if 
a = c*C3, we call the set {1,3} = supp a the support of a. Let X G Vn 
be a partition whose renewal points separate a and b, i.e., there is a 
r G Z such that 

r G 71{X), supp a C {. . . , r — 2, r — 1}, supp 6 C {r,r -\- 1, . . .}. 
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Then 



• Translational co variance : Shifting a partition is equivalent to shifting 
an observable: for every observable a ^ A, 

ujp+x{a) = uJxiT^ia)). 

Here for example, t!^{c2cI) = CgCg. More generally, r^^ is the automor- 
phism on the observable algebra induced by the magnetic translation 
t{i'^/Re^) in the one-particle Hilbert space. Note that the shift trans- 
forms a partition of {0, . . . ,pN — p} into a partition of {p, . . . ,pN}. 

• Locality: let a be a local observable. Let N & N and X, X' G Vjq 
be two partitions. Suppose that X and X' coincide on some interval 
containing supp a, i.e., they have common renewal points r, s such that 
that supp a C {r, . . . , s — 1} and 

11{X) n {r, . . . , s} = 7^(A") n {r, . . . , s}. 

Then 

These properties are, again, a consequence of the product rule Eq. fITBl) . 

For later purpose let Voo be the set of partitions of Z into sets of the 
form {pj,pj + 1, . . . ,pj + pn — 1}, n & N, j G Z. Thus we explicitly forbid 
partitions that contain an infinite component, and every lattice site has a 
renewal point to its left and a renewal point to its right. For X G Voo, and 
a a local observable, let r and s be two renewal points of X enclosing the 
support of a. We can restrict X to a, partition Xr^s of {r, . . . , s — 1} and set 

ux{a) := ujxr^,{a). 

Locality ensures that this definition does not depend on the precise choice of 
renewal points r, s G Tl{X). 

In this way we obtain a family of linear maps a i— ?■ uxia), indexed by 
partitions of Z, A" G Voo, and defined on local observables a. These maps 
inherit the clustering, locality and translational covariance from their finite 
volume counterparts. They will allow us to write the infinite volume state 
as an integral, 

^^^|^a..(a)dPW 
in terms of a suitable probability measure P on partitions of Z. 
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5 Peierls type argument and consequences 



As explained in the introduction, the algebraic structure described in the 
previous subsections allows for a considerable simplification of the problem of 
thermodynamic limits: instead of looking at a full quantum-mechanical state, 
we can look at a simpler probability measure Pat. A further simplification is 
that Pat can be shown to have a stationary limit P if a certain condition on the 
asymptotics of normalization constants is satisfied, see |JLS] and Sect. 15.41 

The aim of this section is to prove that this condition is indeed satisfied. 
This is shown by viewing the normalization as a partition function for a 
Coulomb gas. In Sect. 15. 2^ we adapt the methods of |AJJj to prove that with 
positive probability, the system splits into neutral subsystems with finite 
interaction, and deduce the required statements on asymptotics (Sect. 15.11) . 

We also derive an auxiliary bounds on correlation functions, needed for 
bosons (even p), in Sect. 15.31 

Throughout this section we fix A = Z, i.e., all integrations are on the 
infinite cylinder, and || ■ || = || ■ refers to the norm on the infinite 
cylinder. 

5.1 Normalization on an infinite cylinder 



As was shown in |JLS] . on thin cylinders (7 large), /x(p, 7) is certainly finite. 
However, what happens on thick cylinders was left open. Here we answer 
this question: /i(p, 7) is finite, no matter how large the cylinder radius is. 

Theorem 5.1. For all p and'-/, the quantity fi{p,'~f) is finite. 

The proof of the theorem uses a form of submultiplicativity, which com- 
plements the supermultiplicativity Cn+ai > CnCm- 




(27) 



exist and are finite, with l//i(p, 7) = not yet excluded: 



< r{p,'j) < 1, < /i(p, 7) < 00. 
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Proposition 5.2 (Submultiplicativity). For all N, M and a suitable constant 

Cjv+m<c(p,7)CjvCm. (28) 
The proof of the submultiphcativity is deferred to the next section. 



Proof of Theorem I5.il All we need to show is that for all and some e > 
0, CArr(p, 7)^ > e. But this follows from Eq. ( 128|) and (127|) . by standard 
arguments: repeated application of Eq. (125]) yields 

C^iv<c(p,7)™-'C]^. 
We take the logarithm, divide by mN, and let m — > 00. This gives 

- log r{p, 7) < ^ log c(p, 7) + ^ log Cn- 

It follows that CNr^Py'j)^ > l/c(p, 7) for all and 7) < c(p, 7) < 
00. □ 



5.2 Particle excess function 

Recall from Eq. (|3]) that I'if . . . , zn)\'^ favors particle abscissas Xk = 
{k — l)p7, k = 1, . . . ,N. Think of the infinite cylinder as a collection of A^ 
annuli of width py, centered around those optimal abscissas, and a left and 
right tail. Let x be at the boundary between two annuli, i.e., 

x = {k- l/2)p7 

for some k = 1, . . . , N — 1. The state should prefer configurations with k 
particles to the left of x, and N — k particles to the right of x. Deviations 
from this optimum configuration are measured by the particle excess function 

K{x; zi,..., zn) ■■= i^{j \xj<x}-k = N- k- \ xj > x}. 

Theorem 5.3. Let P/v be the probability measure on with density pro- 
portional to l^^Ni^i, • • • 5 -2Af)P- Let X and K{x; z) as above. Then, for some 
constant c{p, 7) < 00 independent of N and x 

Pn{K{x;zu...,zn) =0) > l/c{p,^) > 0. 

Prop. I5.2l will follow from a similar lower bound on the probability that not 
only the particle excess vanishes, but moreover particles do not accumulate 
too close to X, see the definition of regular configurations f2reg below. The 
next theorem will be useful for the bosonic wave function. 
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Theorem 5.4. With the notation of Theorem \5.3l For suitable constants 
C,c> and all n eNq 



Pn{\K{x; zi, . . .,zn)\ >n)< Cexp( 



-cn 



3\ 



C and c do not depend on N or x. 

The proofs follow ideas from [A J Jj . They are best understood in the 
light of the plasma analogy: we think of |\l/Arp as the Boltzmann weight for 
a classical Coulomb system. Note that Theorems 15.31 and 15.41 only involve 
the modulus of the wave function and they do not refer to the full quantum 
mechanical setting. 

Proof of Prop. \5.2 . Let N, M E N. Taking into account proportionality con- 
stants, Eq. 03]) becomes 

(27rv/vr/7)-(^+^) 

l^^N+M^zi,. . . ,zn+m)\ = (iv + M)! exp{-U{zi,...,ZN+M)) 

The "energy" U{z) is, up to a multiplicative and an additive constant, the 
energy of some classical quasi ID Coulomb system. It is the sum of a one- 
dimensional energy and a correction, 

U{Zi,...,Zn) ■.= U^{Xi,...,Xn) + V2{Zj-Zk) 

where, for xi < ■ ■ ■ < xat+m, 

N+M 2 
. . . , Xn+m) ■■= {^3 " (-^ ~ ^)^^) 

^2(2) = — 2j9ln|l — exp^— 7(|x| + . 

Let 

n:={ze C^+*^ I xi < ■ ■ ■ < Xn+m, 0<yj< 2ttR} 

be the configuration space for N + M particles on the infinite cylinder, la- 
beled from left to right. We imagine the cylinder split into two half-infinite 
cylinders, separated by x := {N — l/2)p7. Let Qrcg C fl he the set of config- 
urations that are regular in the following sense: 

(i) Each particle is in the "correct" half- cylinder, 

xi < ■ ■ ■ < xn < {N - ^)p7 < xn+1 < ■ ■ ■ < xn+m- 
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(ii) Particles do not accumulate at the boundary: for all j = 1, . . . , N + M, 









)P7 


> 


3-N -- 



2 ■ 



In Jlreg, the system's energy U{zi, . . . , zn) is the sum of energies associated 
with the subsystems of and M particles, plus an interaction term that is 
lower bounded by 

N N+M 

j=l k=N+l 
N N+M 

> -2p'^ ^ ln[l + exp(-7|xj - Xfcl)] 

j=l k=N+l 

>_2p(^l-exp(-p7V2))'', 
for all N and M. As a consequence, 

(27rv^/7)~(^+^^) / exp(- U) < const Cn Cm, 

for some suitable A^, M independent constant. Eq. (1281) follows immediately 
once we know that regular configurations have positive probability, i.e., 

/ exp(-t/) < const / exp(-f/), (29) 

uniformly in A^ and M. 

So it remains to show Eq. fl29|) . To this aim we use a Peierls-type argu- 
ment: with each irregular configuration we associate a regular configuration 
that has a smaller energy. This is done by shifting particles that are in the 
wrong half- cylinder, or too close to the boundary, closer to their optimum 
positions (j — l)p7. Some technicalities arise because the implementing map 
T : — )■ r^i-cg is not one to one and has non-trivial Jacobian, resulting in an 
entropy loss that has to be carefully evaluated. 

We start with the definition of the map T : — fireg. Let 2 G fi. Let 
j < A^. We shift the j-th particle as follows: 

• if Xj > X (irregular particle, wrong half-cylinder), then 
x'j := (j - l)p7 - exp{-\xj - x\). 
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• ii X — [N + 1/2 — j)p7/2 < Xj < x (irregular particle, too close to the 
boundary), then 

• ii Xj < X — {N + 1/2 — j)p7/2 (regular particle), then x'^ := Xj. 

Particles belonging to the right half-cylinder, j > + 1, are shifted in an 
analogous way. To complete the definition of T, let y'j := yj and 

T(z) = (2^(1), . . . , z'^(N+M)) 

with a a permutation that reorders the shifted particles from left to right. 
One can check that T maps Q into Q^cg- 

It is convenient to group together configurations that have the same irreg- 
ular particles. Thus, for J C {1, . . . , N + M} we let Qj be the configurations 
that have x'j 7^ xj, if and only if j G J. This gives a partition of ^2. 

JC{1,...,N+M} 

We claim that for suitable constants ci, C2, independent of A^, M and J, 

[ exp{-U)<c^exp(^-c,J2iJ-N)') [ exp(-f/). (30) 

Eq. (129|) follows from Eq. (l30l) by a summation over subsets J C {1, . . . , + 
M}, noting that 

J2 exp(-C2(j - iV)2) < [11(1 + exp(-C2A;2))'j < 00. 

JC{1,...,N+M} \k=0 J 

We are finally left with the proof of Eq. (1301) . It is shown with the help of a 
change of variables. We refer the reader to |AJJj for the details and content 
ourselves here with some hints about the necessary estimates: 

Entropy: One has to give an upper bound of the maximum number of 
preimages |T"^(z')| of points in z' G T{Qj), and a lower bound for Jacobian 
of |dT(z)/dz| for z & Qj; note that T is almost everywhere different iable. 
The bounds depend on the set of irregular particle labels and are bad when 
there are a lot of irregular particles, but this is compensated by a gain in 
one- dimensional energy AU'^. 
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Energy: The simplest to estimate is the decrease in ID energy, 



AU' := U\z) - U\T{z)) > hY^U - N)' + h 



-X? 



j£j j- wrong 

half-cyl. 



for some suitable constant ki. 

For the V2-interactions, note that only the interactions affecting irregular 
particles are changed. Now, instead of estimating the change in V2-interaction 
directly, we adopt a three step procedure. First, we drop the affected V2- 
interaction altogether, using 

N+M 

Y Y ^2(2,- -^fc) > -hiJi^ 

jeJ k=i 

for some constant k2- Next, we shift Xj h- )■ x'j. Finally, we reinsert the 
dropped interactions with Jensen's equality: for fixed x 7^ 0, V2 is a harmonic 
function of y, hence J V2{x + iy)dy = and 

„ „ N+M 

ldy'< exp(^-Y Y V2{z',-z',))dy'. 

j&j k=i 

All y-integrations go from to 2'kR. At the end, we obtain an estimate on 
^/-averaged Boltzmann weights, 

j exp(-f/(2))d2/<exp(A:2|Jp-A?7i) j exp(- f/(T(z)))dt/', 

with J dy' the integration over ^/-coordinates of T{z). Eq. fl30l) is obtained 
by combining this last estimate with the entropy estimates. □ 

Proof of Theorem \5.3[ In the notation of the proof of Prop. 15. 2[ regular con- 
figurations have zero particle imbalance and, therefore, Eq. f l29|) gives 

Pn+m{K{x; Zi,..., Zn+m)) > PN+ui^reg) > V coust > 0, 

uniformly in and M. □ 



Proof of Theorem \5.4\ Again, it is enough to have a closer look at the proof 
of Prop. 15.21 We remark that if the particle excess is negative, there are 
too many particles in the right half-cylinder: if K{x; Zi, . . . , zm+n) < —k 
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for /c G N, then > ■ ■ ■ > x^-k+i > x. Thus the particles with labels 
— A; + 1, . . . , are irregular. Since 

fc-l / oo \ 

E exp(-C2^(j-iV)2) <exp(-C2 5^.2) 11(1 + ^^P(-^2f )) 
jc{i,...,]v+A/} jeJ s=o \i=o / 

{JV-fc+l,...,JV}C J 

we find that the weight of configurations with particle excess < —k is of order 
exp(— const /c"^), uniformly in and M. A similar reasoning can be applied 
to positive particle excess. This proves Theorem 15.41 □ 



5.3 Uniform moment bounds for lattice bosons 

Before we investigate the thermodynamic limit, we give bounds on correlation 
functions in terms of the ID lattice system. This is an auxiliary result needed 
for bosons {p even) only. For fermions, the bounds given below are trivial 
since there can be at most one fermion per lattice site. 

Recall that c^,Cfc, G Z, are the creation and annihilation operators for 
the orbital ipk{z:). Also, we denote {a)N '■= {'^N,ci'^N)/CN, with integration 
on the infinite cylinder. 

Proposition 5.5. Let n E N. Then for some suitable constant D > 0, all 
n,N E N and all ki, . . . , kn, mi, . . . , m„ G Z, 

\{nk,---nkjN\<D'' (31) 

\{cl---clCm,^...Cm,)N\<D\ (32) 

Proof. We start with the proof of the first inequality in the case ki = ■ ■ ■ = 
kn = k. Thus we seek to estimate {n'^)N- Let A := [a,b) x [0,2-7ri?] be the 
annulus of width 7, centered at x = k'j, i.e., a,b = k'j ± 7/2. For a given 
configuration (21, ... , z^), the number of particles in A is 

Na{zi, ...,zn) = 1 + K{b; zi,...,zn)- K{a; zi,...,zn). (33) 

It follows that Na has moments of all orders, which by Theorem 15.41 are 
uniformly bounded, E^{N2) < w"" for all n,N,A and some suitable w > 0. 
Here refers to the measure Pjy on Z'^ with density oc l^^^izi, . . . , z^)]"^, 
and we use that k^ exp{—k^) grows at most exponentially in n. In order 
to translate bounds from the continuum picture to the lattice picture, we use 
the formula for the n-point correlation function at Zj = z'j-. 

Pni^ly ■ ■ ■ y ^n] Zi, . . . , Zn) = ^^(Cfc^ " " " Cfc^Cm„ " " " Cmi) N 

ijjkiiZl) ■ ■ ■ i^kr,{Zn) i^miiZl) m„ [Zn ). (34) 
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The sum is over the integers kj^rrij G Z. When we integrate over the y- 
coordinates, the contributions with kj ^ irtj vanish and the non-vanishing 
contributions are positive. Hence, integrating over A", we obtain 



E 



M [Na{Na - 1) ■ ■ ■ (iVA - n + 1)) > n(7)"((4)"(c,)"> 



TV 



The constant ^(7) < 1 comes from integrating the Gaussian over the 

annulus A; it does not depend on or k. For the left-hand side, we have 
used that the diagonal n-point correlation is a factorial moment density (see 
|DVJt Sect. 5.4] for an explanation of this notion). The right-hand side is a 
factorial moment too, since 

{cinckT = M^k - 1) ■ ■ ■ (rifc - n + 1) =: . 

Now, moments are linear combinations of factorial moments with positive 
coefficients Aj „ > 0, the Stirling numbers of the second kind |DVJt Sect. 
5.2]. We deduce 

n n 

i=i i=i 

It follows that 

for all n. Repeated use of the Cauchy-Schwarz inequality | (a*6)| < {{a*ay^'^{b*by^'^ 
yields Eq. f l3T|l . Since factorial moments are smaller than moments, it follows 
that Eq. fl32l) holds when ki = mi, . . . ,kn = rrin, from which the general case 
kj 7^ rrij is deduced, again, with the help of Cauchy-Schwarz. □ 

As an application we prove Prop. 12.41 



Proof of Prop. \2.4\ when A = Z. We know that {cl.Cn) at 7^ whenever k ^ n, 
see |JLSj . thus 

pN—p 



k=0 

For every integer k, we can factorize 



\iJk{z)tpk{z')\ oc exp(^-i((x - k-ff + {x' - k-ff)^ 

= exp^— -(x — x')^^ exp^— -(x + x' — 2k'-fY^. 

By Prop. 15.51 the occupation numbers are uniformly bounded. The proof is 
concluded by noting that Xlsez ^^P(~(^ ~ ■^7)^/'^) bounded too, uniformly 
in a; e M. □ 
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5.4 Stationary renewal process 

Recall the weights p7v(A') and the corresponding measure Pat on Vn, see 
Eq. (124|) . When we integrate on the infinite cylinder, the translational co- 
variance of the functions Vx{zi, . . . , Zi^(^x)) gives Eq. (14. 3 p the form 



Cm = \\^n\?= Y1 



X • • • X Q^^^. 



niH \-no=N 

As shown in [JLSj . the positive numbers «„, relate to the quantities 7) 
and /i = 7) from Section EH] through 

00 00 

^ a„r(p, 7)'' = 1, na„r(p, 7)" = 7). 

n=l n=l 

Thus we may consider pn '■= a„r(p, 7)*^ as a probability distribution on N, 
with finite expectation 7). We have 

Mat := CArr(p,7)^ = ^ PniX---Xp„^. 

niH ynB=N 

This is the probability for a renewal process with waiting time distribution 
(pn) to have a renewal point at A^, given that it had one at 0. For our purpose 
it is more convenient to view (pn) as a measure on pN. The renewal process 
is then called p-periodic. 

The weights p]\f{X) from Eq. become 

Pn{X) =Pni ■■■Pnu/UN 

where X is & partition of {0, 1, . . . ^pN — 1} into D consecutive intervals of 
lengths prii, . . . ^pno. Thus Pat is a p-periodic renewal process conditioned 
on and pN being renewal points. 

Now, given a distribution on N with finite expectation /i, there is a stan- 
dard way of defining a stationary renewal process, or, in our case, a p-periodic 
renewal process. To each partition of Z we associate the indicator function of 
the renewal points. In this way Poo inherits the product topology and Borel 
cT-algebra from {0, 1}^. There is a unique measure P on Voo such that for all 
integer j, 

P(p?' is a renewal point) = — r, 

and for all j < k, 

is a renewal point | pj is a renewal point) = Uk-j- 
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This measure P is invariant with respect to shifts by multiples of p. The 
following three lemmas will be fundamental for the investigation of the ther- 
modynamic limit. 

Lemma 5.6 (Pat — ?■ P). Let r, s G N have distance > d from and N, 
d < r, s < N — d. Consider the event Srs that pr and ps are renewal points of 
A" and moreover the restriction of X to {pr, . . . , ps — 1} coincides with some 
given partition of this subset. Then 

|PAr(^^rs) - P(^^rs)|< COUSt PAr(£:,.^) SUp \Uk - fJ^~'^\. 

k>d 

Proof. Suppose that we require that X coincides, in {pr, . . . , ps — 1} , with a 
partition into successive intervals of lengths pni, . . . , priD, ni+- ■ ■+ni:, = s—r. 
Then 

P(£„) = ■■■Pno 
^N{£rs) = UrPm " " -PnoUN-s/uN 
Pjv(g,,) - ¥{Srs) ^ ^ UN 

FNiSrs) HUrUN-s 

The claim then follows from the observation that m„ — t- fi~^ as n — )■ oo. □ 

Lemma 5.7 (Long intervals are unlikely). Let a, /3 G Z with (5 — a>pd for 
some d G N and c{p, 7) > as in Prop. 15.21 Then, for all G N, 

Pat (A" has no renewal point in {a, . . . , /3 — 1}) < c{p, 7) kpk- 

k>d 

Proof. The submultiplicativity Prop. \572\ gives 

— — < — < c{p,i). 

UN Un 

The probability that a partition X contains the interval {pj, . . . ,pj +pn — l} 
is therefore 

UjPnUN~j~n/uN < c{p,-f)pn. 

We have to sum over pairs (j, n) such that pj < a and pj + pn > /3. Noting 
that 

X] Pfe = - d)pn < ^npn, 

j>0 n>j+d n>d n>d 

we obtain the desired inequality. □ 
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Lemma 5.8 (The renewal process is clustering). Let < r < s < have 
mutual distance s — r > d. Let C be the event that pr is a renewal point and 
moreover the restriction of X to {0, . . . ,pr — 1} coincides with some given 
partition of this subset. Define 7^ in a similar way referring to the subset 
{ps, . . . ,pN — 1}. Then 

\Fn{C n 7^) - FNiC)FN{n)\< const P7v(>C)P7v(7^) sup \uk - fi~^\. 

k>d 

Proof. The proof is similar to the proof of Lemma 15. 6^ except that the term 
to be estimated, in the end, is 

{un-s/un) {Ut/un) UN-sUr 

Again, this term is small for large d because u„ — fi~^ as n — )■ oo. □ 



6 Correlation functions 

In order to prove our main results, all there is left to do now is to go from 
the convergence of measures Pjv — )■ P to convergence of states {a)^ — )■ (a). 
This is essentially an interchange in the order of summation and limits, and, 
therefore, involves some technical estimates; we hope, however, that we have 
conveyed the simplicity of the underlying idea. 

The proofs follow ideas from [JLS] , even though the presentation is slightly 
different as no use was made, in that work, of the framework of quasi-state de- 
compositions. Some estimates become more involved because we extend the 
results to bosons, so that creation and annihilation operators are unbounded. 
Another novelty is the proof of the insensitivity towards the precise choice 
of the domain of integration in Sect. 16.41 

6.1 A variant of the Cauchy-Schwarz inequality 

Recall Eq. ([25]) 

{a)N = -^{'^N,a^N) = V] PN{X)uJx{a). (35) 

We would like to deduce from Lemma 15.61 that 

(a)^^ [ d¥{X)uM=-{(^)- (36) 
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To this aim we will need to control the contribution of unfavorable partitions 
X. When a is a diagonal, bounded operator, and J-" is some collection of 
partitions, we can write 

J2PNiX)Ma)\<\\a\\¥NiX eJ'). (37) 

For observables that are not diagonal or bounded, we use a variant of the 
Cauchy-Schwarz inequality Tr(pA*5) < {Ti pA*Ay/^{Ti pB*By/^. This will 
help us, in some situations, replace the estimate (1371) by an estimate of the 
type 

for some a-dependent constant d[a). 

Lemma 6.1. Let p be a trace class, positive operator in l?{H) with ma- 
trix {pij)ij^s. Let A = (Ay)jjgN be an operator in £^(N). Suppose that 
Trpyl*yl < oo, with A possibly unbounded, and that A has at most d{A) 
non-zero matrix elements Aij per row. Then for all S* C N, 



im jes 



Proof. It is enough to treat the case when A has at most one non-zero matrix 
element per row. Write pij = J2neN ^nX^x^ with Ai > A2 > ■ ■ ■ > the 
eigenvalues of p and G i^(N) the normalized eigenvectors. There is a map 
: N — N assigning to each row i the column j = with the non-zero 
entry. If all entries in the row vanish, is arbitrary. Then 

E E -EE ^«l^i^J:<^(i)^^(i)l 

<E^»(Ei-'.f)'"(I^i-*'*)=^*(')if" 
< E^«Ei^"n E^"EiA0w^;»r 

\ n ieS ) \ n i^S 

\l/2/ \l/2 

J^Pii) (TrpA^A) . D 



< 



As a first application of the previous lemma we show: 
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Proposition 6.2. Let a = cl ■ ■ ■ c\^Cq^ ■ ■ ■ Cg^ be a Wick ordered monomial 
of creation and annihilation operators. Then, for all N eN, 



J2 Pn{X)\coM\ < {o.*a) 



Proof. For every m = (mi, . . . , tun) with < mi < • • • < tun < pN — p, 
there is at most one m' such that {ipm', CLipm) 7^ 0. Thus we are in the setting 
of case 2. of the previous lemma, with d{a) = 1, and we deduce 



As a consequence, the candidate hmit (a) in Eq. (1361) is well-defined: 

Corollary 6.3. Let a be a linear combination of products of creation and 
annihilation operators cl, Cm- The integral from Eq. is absolutely con- 
vergent. 



This is shown by passing to the limit — )■ 00 in the previous proposition, 
and using Lemma [5.61 and the uniform moment bounds from Prop. 15.51 

Remark. Cor. 16.31 has an analogue when a is a local, bounded operator, for 
example a = exp(ic^ci). This can be shown with the help of yet another 
variant of the Cauchy-Schwarz inequality. 

Note that despite the notation (■), we do not know yet whether the linear 
map a (a) is positive or defines a proper state. 

6.2 Thermodynamic limit 

In this section we prove Theorem 12.11 when the infinite cylinder A = Z is 
chosen as the domain of integration. We first estimate the difference between 
left- and right-hand side in Eq. fl5Bl) . 

Lemma 6.4. Let a be a Wick ordered monomial supported in {a, . . . , /3 — 
1} C Z. Let d E N and Sd be the set of partitions with renewal points in 
both [a—pd,a] and [/3, /3 + pd]. Then 



J2 Pn{X)\uM\ < 



1 




^ PN{X)\uJx{a)\ < (a*a)!v/^ (^2c(p, 7) ^ /cp^) 



1/2 
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Proof. We abbreviate m = (mi, . . . , mjy) and 



aAr(m) = aN{mi, . . .,171^), ij{m)=tprm{zi) ■ ■ ■ V^m^r(^Af)• 
Let p{m,m') := aM{m)aN{m')/CN- For d G N, let S* be the set of m's 
without renewal points in [a — pd, a] or without renewal point in /3 + 
For every m, {■ip{m'),atlj{m)) is non-vanishing for at most one vector m'. 
Thus by Lemma f6.lt 

PN{X)\u;x{a)\ < j;^^ ^aN{m)aN{m'){tlj{m),aij{m')) 
< {a*a)f [j2Mm)\ViN\a ^^"^ 



■'N 



and we conclude with Lemma 15.71 □ 

Lemma 6.5. Let a be a Wick ordered monomial supported in {a, . . . , /3 — 
1} C Z, (i e N, and Sd as in the previous lemma. We may consider Sd as a 
subset of Vn or of Voo- Then 



< const ^sup — /i pj^f{X)\ujx{(i) 

X&Vn 



-k>d 



This is a consequence of Lemma [5.61 and the locality of the ux- 

Proof. A partition X & Sd has a renewal point r G [a — pd, a] and a renewal 
point s G [f3, (3+pd]. Choose r the largest possible and s the smallest possible. 
The value of uJxia) depends only on Xr^s, the restriction of the partition to 
the volume enclosed by the renewal points. The contribution of partitions 
with the same r, s and same restriction Xj.^s is of the form ux^ ^i^WNi^rs)- 
A similar form can be derived for partitions of the infinite lattice Z. By 
Lemma 15.61 



UJx.Ja)FN{Srs) - iOx^Ja)F{Srs 

< constfsup \Uk - fi'^l) \^^Xr^s{(^)\ ^N^Sr 



^k>d 

and the proof is concluded by summing over r, s and Xrs- □ 



36 



The previous two lemmas, together with Prop. 16.21 and Lemma [5.61 yield 
a lattice version of Theorem 12.11 

Corollary 6.6. Let a be a Wick-ordered monomial whose support has dis- 
tance > pd to the lattice boundaries, supp a C {pd, . . . ,pN — pd — 1}. Then 

1 /2 

{a)N-{a) < const (sup (a*a)M) ( 4c(p, 7) fcpfe)^^^ + sup - ) . 

^ M ^ ^ 7-^ k>d ^ 

k>d — 

Note that the upper bound goes to as c? — 00 because ^ kpk < 00, 
Uk and because of the uniform moment bounds from Prop. 15.51 

Proof of Theorem \2.1\ when K = Z. In view of Cor. 16. 6[ all there is left to 
do is to go from lattice correlations to continuum correlations. This is eas- 
ily achieved, thanks to the explicit relation between lattice and continuum 
and the Gaussian decay of the one-particle functions ijjk{z)- For the sake of 
clarity we write down the proof for the two-point functions only; the other 
correlation functions can be treated in a similar way. Recall 

p^'^ {Zi, Z2; Z[, = {c*^C*^CeCk)Nlpk{zi)iJe{z2)llJm{z[)lpn{z'2). 
0<k,£,m,n<pN—p 

(38) 

Therefore we define 

p2{Zl,Z2;z[,Z2) := {c*mCnC£Ck) ^k{Zl) M^2) iJm{z[) ^n{z'2) 

k,i,m,n& 

with the bulk expectation (■) as in Eq. fl36l) . The infinite sum converges be- 
cause the expectations (■ ■ •) appearing in the sum are bounded as in Prop. [531 
and 

sup iV^ml^)! oc sup exp[— (x — < 00. 

The translational covariance of uJx{o,) and the stationarity of the renewal 
process P make the two-point correlation function pi"^/ i?-periodic in the sense 
of Eq. dS]). 

The difference between the finite volume and infinite volume two-point 
function is a sum over integers k,l,m,n. Suppose that zi, z[, Z2, Z2 are at 
distance > D > 2pd'-f from the boundaries of the cylinder x = and x = 
pN'y. The contribution to the two-point function from summands with all 
of the four indices between d and pN — pd can be bounded with the help 
of Cor. 16.61 The contribution from quadruplets [k, I, m, n) with k < pd, is 
bounded by a constant times 

Y \M^'i)\ oc J2 exp(-K - k^Y/2) < 5^exp(-p2(rf + «:)V/2). (39) 

k<d k<pd K>0 
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which is small for large d. Contributions where another index is smaller than 
pd or larger than pN — pd can be bounded in a similar way. □ 

6.3 Clustering 

The clustering for A = 2 is deduced from the renewal process Lemma 15.81 
in the same way as the thermodynamic limit was deduced, in the previous 
section, from Lemma 15.61 Again, we start with lattice correlations. 

Consider a and b two Wick-ordered monomials of creation and annihi- 
lation operators with supports at mutual distance > 3pd for some d E N. 
Thus let G Z such that P — a > 3d and suppa C {. . . ,pa — 1} and 
supp b C {p(3, 0/3 + 1, . . .}. Let J-" be the set of partitions X with a renewal 
point in \pa,pa + pd] and Q the set of partitions with a renewal point in 
\p/3 - pd,p(3]. 

We write the expectation of ab as a sum over partitions in HQ plus a 
remainder; similarly for a and b. In order to estimate {ab)^ — {o)N{b)N, we 
need to estimate four terms, the three remainders, and the difference 

J2 PNiXPM^xib) - (j2pM{X)uj;,{a)') (j2PN{X)ujx{b)y (40) 

Recall that if X E (IQ, then uxiab) = ux{(i)u}x{b)- The remainders for a 
and b are simplest to estimate, 

1 /2 

^PN{X)\uJx{a)\ < (^{a*a) Nc{p,-f)^kpk^ , 

X^T k>d 

and an analogous inequality holds for b and Q. The proof is similar to the 
proof of Lemma 16.41 The difference (140|) is bounded by some constant times 

(sup |/i"-^Mfc - l\){a*a)]^'^ {b*b)]^'^. 

This is shown with the help of Lemma 15.81 and proceeding in a way similar 
to the proof of Lemma 16.51 

The remaining estimate of the contribution to {a\})N from partitions not 
in J-" n ^ is slightly more involved. First we switch to an occupation number 
picture. With m = (mi, . . . , mAr), we associate a sequence of occupation 
numbers n = (no, rii, . . . , UpN^p) in the obvious way; for example, when N = 
2, p = 2, and (1711,1712) = (1, 1), we have (rio, ni,n2, n^) = (0, 2, 0, 0). We let 
|n) be the normalized wave function proportional to the (anti-)symmetrized 
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product of iprrni----, i^ruN- The |n)'s form an orthonormal system. The many- 
particle wave function becomes 

^Af = ^AN{no,ni, . . . ,npN_p)\noni ■ ■■UpN-p). 

n 

When p is odd, the coefficients AM{n) are in one-to-one correspondence with 
the coefficients ajq{m). When p is even, the correspondence is up to factors 

Next, we observe that the notion of renewal point, since it does not depend 
on the order of the mj's, can be transferred to occupation numbers. Hence 
pk is a renewal point of n if and only if 

pk—l pk~l 

'^Uj = k and ^ jnj = pk{k - l)/2. (41) 

j=0 j=l 

The admissibility condition for m leads to the following property, valid when- 
ever A7v(n) 7^ 0: for all A; = 1, ■ ■ ■ , - 1, 

pk—l pk~l 
j=0 j=l 

Now let S be the set of n's with no renewal point in [a, a + pd] or no renewal 
point in [/3 — pd, f3], and E the set of pairs (n, n') with no common renewal 
point in [a, a -|- pd\ or no common renewal point in [f3 — pd, f3]. Thus n E S 
is a sufficient, but not necessary, condition for (n, n') to be in E. Write 

V pN{X)ux{ab) = ^ V AN{n)AN{n'){n\ab\n'). (43) 

The next lemma paves the way for an application of a Cauchy-Schwarz- 
inequality to Eq. fH3|) . 

Lemma 6.7. Let a, b, E as described above, and suppose in addition that 
a preserves the total particle number, i.e., it is the product of n creation 
operators and the same number n of annihilation operators. If (n, n') G E 
and 

AN{n)AN{n'){n\ab\n') 7^ 0, (44) 

then n E S OT n' E S. 
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Proof. Recall the notion of support of an observable: when a = c^^Cg, supp a = 
{1,3}. If Eq. fl44l) holds, then n and n' must coincide outside supp(a) U 
supp(fe). We may assume, without loss of generality, that 

pa— I pa— I 

5] J'^j- > 5Z J^i- (45) 
j=i j=i 

(Otherwise swap n and n'.) If a preserves the total particle number, e.g., 
a = c\c2, then n and n' must have the same number of particles in supp a. 
It follows that for every site v between supp a and supp 6, they have the same 
number of particles to the left of f , 

riQ + ni^ hnt,_i=noH hn^.^, pa<v<p[5-l. (46) 

Suppose that n has a renewal point pk between supp a and supp b, a < k < p. 
Then n has k particles to the left of w = pk, and by Eq. (146|) . so has n'j. 
Moreover, since n and n' coincide outside supp(a6), Eq. (145|) implies that 

pk—l pk—l 

jn'j < jrij = pk{k - l)/2. 
j=i j=i 

Eqs. (1121) and fHTl) then imply that pk is a renewal point of n' too. Thus 
every renewal point of n between supp a and supp b is in fact a common 
renewal point of n and n'. Therefore, if (n, n') G E, necessarily n E S. □ 

Now we can apply Lemma [6.11 which together with Lemma [5. 71 yields 

1 /2 

I J2 PN{X)u;^{ab)\<2{{abyab)]i^(2c{p,^)J2kpk) ■ 

X(^T<r\Q k>d 

Using our uniform moment bounds, we obtain: 

Proposition 6.8. Let a be a product of n creation operators andn annihila- 
tion operators, and b a product of m creation and m annihilation operators. 
Suppose that supp a C {. . . ,pa — 1}, supp 6 C {pf3, . . .} and (3 — a > 3d. 
Then, for some suitable constant Km+n o,nd for all N, 

ab)N - {a)N{b)N < ^m+n( (V" ^i^Pfe)^^^ + sup \Uk - fi^^ 

k>d — 
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Remark (Off-diagonal decay). Suppose, for example, a = cl and b = c^CmCn 
for some k,l,m,n G N. Then {a)iy = = {b)^ = and, because of y- 
momentum conservation, (ab) at = unless k + I = m + n. Imagine shifting 
b along the x axis, 

b c*i^^Cm+dCn+d = r^{b). 
Then {aT^ib))^ ^ unless /c + (/ + d) = m + + 2(i, so we find 

d> k + l-m-n {aT^{b))N = = (a) N{T!^{b)) n- 

The same reasoning applies to higher correlations where a, or b, does not con- 
serve particle number. Hence there is clustering for off-diagonal correlations 
too. 

Proof of Theore'm \2.3[ Just as in the proof of Theorem 12.11 on p. [371 all we 
need to do is pass from lattice (Prop. 16.81) to continuum. Again, for sim- 
plicity we write down the proof only for the two-point correlation. Let Zi, z[ 
have large distance from 22,^2 along the cylinder axis, i.e., xi,x'i < wy and 
2^2,3^2 ^ ("W^ + 3-D)7 for some integer w and large D. We write the sum in 
Eq. (!38|) as a main contribution M, plus a remainder. The main contribution 
consists of those summands where k,m < w + D and l,n > w + 2D. Because 
of Prop. 16.81 we have a bound 

\{c*mCnClCk)N ' {c*mCk) N {clCl) n\ < COnst g{D) 

for some function g{D) — t- as -D — t- 00. It follows that the difference between 
the main contribution M and 

( i^mCk) nM^i) [ Yl «C/)iV^i(2;2)^n(4)) 

k,m<w+D l,n>w+2D 

can be bounded by a constant times g{D) too. But in this last term we 
recognize the main contribution to {zi, z[)p^ {Z2; z'2); so we are left with 
the remainders to estimate. This is easily achieved: the remainders, for the 
one-particle matrix as well as for the two-point function, can be estimated 
by terms of the type fl39l) . □ 

6.4 Non-influence of the domain of integration 

In this section we show that the precise choice of domain of integration 
does not affect bulk correlations. The idea is to rewrite integrals over A 
as integrals over the infinite cylinder Z with the indicator function of A in 
the integrand, and then translate the indicator into a lattice operator. This 
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will allow us to view the indicator function as a quantity that lives at the 
cylinder's boundaries, and to decouple this boundary perturbation from bulk 
correlations with the help of the state's clustering. We start with a simple 
computation. Let mi,...,mjv € Z (not necessarily ordered or distinct). 
Then 

/ \lprm{Zl)\'^ ■ ■ ■ \lpmMizN)\'^dZi ■ ■ ■ dz^ = H^milli" " ' \\'ipmN\\l■ 
We can rewrite this as 

using the diagonal operator in L?'{Z^) 

■Jn,A ■ 'tpm.i ® ■ ■ ■ IpmN ^(11^ 

(we set Jn,a equal to in the orthogonal complement of the lowest Landau 
level, i.e., the space spanned by the ipk{z)). The Fock space version of this 
operator, again denoted Jn,a, is 

pN—p 

|2 \"fc 



jn,a = n (ii^^iia)' 



fc=0 



In this product, only boundary terms, k small or close to pN — p, contribute. 
Indeed, ioi < k < pN - p, and A = [a, b] x [0, 2ttR], 



2 

fcllA 



1 ^-i--k^?d^ = 1 - e, - <5p^_p_fc. (47) 

J a 



The error terms and 6j depend on the precise choice of the domain of 
integration. They are small when k oo, resp. j — ?■ oo. For example, when 

a = 0, b = {pN — p)7, 

1 

Sk = ek = —1= / e'^^^ds {k oo). 

Lemma 6.9. The normalization and the one-particle density for \&Ar with 
domain of integration A are given by 

II^Tvlli = (*^, ^7V,A*7V)Z 
N,Af n {'^N,clJN,ACk'^N)z , / x , / .x / an /aq\ 

Pi {z]z)= > — -— Mz)Mz'), {z,z gA). (48) 

^ \^ N, 'Jn,a^ n)z 
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Similar formulas hold for ra-point correlations. 

Proof. The formula for the normalization is a consequence of the computa- 
tions used to define Jn,a- For the one-particle matrix, we note 

= {N-iyr^^^ ^ aN{k,m2, . . . ,mN)iljm^{z2) ■ ■ -ipmNizN)- 

m2,...,mM 

Let d{m2, rriN) := \ \^.m2 1 Ia " " " Wi'mN I Ia- We obtain 

(^AT, 4 JN,KCk^N)z = \aN{k,m2, . . . ,mN)\^ d{m2, . . . ,mN)- 

m2,.--,"iJV 

On the other hand, times the integral of "^n^z, Z2, . . . , zn)'^n{z', Z2, ■ ■ ■ , z^) 
with the Zj integrated over A, equals 

iV-^ ^ \a'N{mi, . . . ,mN)\^d{m2, . . . ,mN)ipmi{z)iprni{z'), 

mi 

and the proof is easily concluded. □ 

Because of Eq. f H7|l . it is natural to think of our lattice indicator Jn,a as 
a product of a left, bulk, and right term. We write Jat^a '■= CBTZ with 

d-l 

£ = ]^(l - - 5pN~p~j) \ 

B a similar product for j from d to pN — p — d, and TZ the product from 

pN — p — d to pN — p. 

Lemma 6.10. For all and d with pN — p > 3d and a suitable, A^- 
independent function f{d) with /(d) — )■ as d — ?■ oo: 

sup I (4 JAr,ACfc)7V,2 - {C)n,Z (clCk) N,z{T^) N,z\ < f{d), 
d<k<pN-p-d 

\{Jn ,a)n,z — {^)n,z {TZ)n,z\ < fid). 

Proof. Let d be large enough so that < 1/4 and 5^ < 1 /4. Let c > such 
that ln(l — x) > —cx when |x| < 1/2. For pN — p > 3d, we have 

pN—p—d 

1 > B > 1 - c ^ (ej + SpN-p-j)nj, 
j=d 
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as an operator inequality. Noting that B and c^, commute or anticommute 
with C and TZ, 

\{ciCBnck)N,z - {clcnck)N,z\ 

I \ 1/2 / \ 1/2 

< (£4(1 - By^L) (7^4(l - BYkH) 

\ I N,Z \ I N,Z 

2 1/2 pN~p-d 2 

<((4(l-i3)c,) )^^<((c 5^ (6,+V-,_,)%) )^^. 

Using the uniform moment bounds from Prop. 15. 5^ this can be further upper 
bounded by a constant times YlJLdi^j + ^i)' which is finite and goes to as 
d — )■ oo. 

Next, we want to use Prop. 16.81 in order to decouple C, c\ck = and 
B. We note, first, that Prop. 16.81 is not directly apphcable to £ and B, since 
these operators are not polynomials of creation and annihilation operators. 
Power series expansions yield a simple remedy. To illustrate the procedure, 
we explain how Prop. 16.81 can be applied to a product of two exponentials. 
We have, for suitable K,D > and e{d) — )■ as d — > oo. 



(exp(-srio) exp(-trid))^ - (exp(-sno))^(exp(-tnd)) 
k\ q 



N 



< K^J^T^-J^''^''^^^) ^ Kexp{sD)exp{tD)e{d). 



k=0 q=0 

Here we have used Prop. 15.51 in the form {n^fi^ )T < 0^+" for suitable D. 
Something similar can be done to bound 

{Chkn)N,z - {jC)N,z{f^k)N,z{'^)N,z. (49) 
The upper bound will involve 



JpN-p-j, 

j=0 

and a similar term for the right boundary. We note that this term can be 
bounded, uniformly in and d, and deduce that the absolute value of ( 149|) 
is bounded by some function g{d) with 5f((i) — )■ as c? — )■ oo. In combination 
with our earlier bound which justified the replacement i3 ~ 1, this proves 
the first inequality of the lemma when d > d^ for some do. When d < d^, we 
note that the left-hand side of the inequality can be bounded, uniformly in 
N, and set f{d) equal to that bound. 

The proof of the second inequality is similar. □ 
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Next, we observe that {C)n^z and {TZ)n,z stay bounded away from as 
— 7- oo. Indeed, if (tfc) is a sequence of numbers in [0, 1] and Ylk^k < oo, 
Jensen's inequality gives 




j=0 j=0 



oo 



> expf/s:^ln(l > 0. 



j=0 



Here i^' is a uniform upper bound for the occupation numbers {nk)N,z- This 
argument can be adapted without problems to lower bound the expectations 
of C and TZ. Note that, as operators with norm < 1, they have expectations 
upper bounded by 1. 

As a consequence, we can pass to quotients and deduce from Lemma 16. 101 



for some A^-independent function g{d) which goes to oo as d — ?■ oo. From here 
the proof of Theorem 12.11 for the one-particle matrix for general A, i.e., the 
insensitivity to the domain of integration, is proven with the help of Eq. psj) 
by imitating the proof of Theorem 12.11 for A = Z on p. [371 The proofs for 
general n-point functions are similar. 

6.5 Symmetry breaking 

We conclude the paper with a proof of Theorem 12.21 which is essentially a 
consequence of results of [A J J] . Let us also recall that on thin cylinders, the 
slightly stronger statement that the one-particle density (and not just any 
correlation function) has a non-trivial period was proven in |JLSj . 

Proof of Theorem \2.2[ The diagonal infinite volume correlation functions 

) are the correlation functions (= factorial moment 
densities) of some point process on Z. Because of Theorem 12. 1^ the corre- 
sponding measure P is the limit, in a suitable sense and up to shifts, of the 
measure Pn with density oc I^^n^zi, . . . , z^)]"^, choosing the finite cylinder 
—py/2 < X < p{N — 1/2)7 a.s the domain of integration. With this choice 
P/v is exactly the Gibbs measure for particles moving in a neutralizing 
background, studied in [AJJj . 

Therefore, by jAJJl Theorem 3.1], if we shift P by 6* e M along the x- 
axis, we obtain a measure P^ which is singular to P unless 6 is an integer 



sup 

d<k<pN-p-d 



{clJN,ACk)N,Z 

{Jn,a)n,z 



{clck)N,z < g{d) 
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multiple of p7. Now, Theorem 15.41 together with Eq. f l5^ shows that the 
point process satisfies conditions which ensure that it is uniquely determined 
by its correlation functions |DV J] . Remember that passing from correlation 
functions to the point process is like passing from moments of a probability 
measure to the measure itself. Thus if 9 is not an integer multiple of p7, 
the ^-shifted measure must have some correlation function which is different 
from the one for the original measure P. This proves Theorem 12.21 □ 

Remark. Repeated shifts of the infinite-cylinder state 0Ji{-) = (■) by i'^/R in 
the x-direction yield states U2, ■ ■ ■ ,ujp. By Theorem 12.21 those p states are 
distinct. They are actually also disjoint (this notion generalizes mutual sin- 
gularity of probability measures). This follows from general arguments [BRl 
Sect. 4], combining the fact that cui, . . . , are distinct and mixing (by The- 
orem [23]), hence ergodic with respect to shifts in the x-direction. 
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